Conductance of telescoped double-walled nanotubes from perturbation calculations 
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In a telescoped double- walled nanotube (TDWNT), with the inner tube partially extracted from 
the outer tube, the total current is forced to flow between the layers. Considering the inter layer 
Hamiltonian as a perturbation, we can obtain an analytic formula for the inter layer conductance. 
The accuracy of the perturbation formula is systematically improved by including higher-order 
terms. The inter layer interaction effective in the perturbation formula is the product of the inter layer 
Hamiltonian and the wave function. It clarifies the effects of the spatial range of the inter layer 
Hamiltonian and the band energy shift. 



I. INTRODUCTION 

Electronic and mechanical properties of single-walled 
nanotubes (SWNTs) originate from the tt and a bonds 
in their honeycomb lattice i"— Although the inter layer 
bonds are considerably weaker than the intralayer tt and 
a bonds, they are important for the formation of double- 
walled nanotubes (DWNT), multi- walled nanotubes, and 
nanotube bundles. The telescoped double-walled nan- 
otube (TDWNT) shown in Fig. 1 was formed from a 
DWNT by partially extracting the inner tube. Since the 
rigid honeycomb lattice is relatively unaffected by the 
weak inter layer bonds, the basic motions of the TDWNT 
are limited to the interlayer motions caused by relative 
slide and rotation between the outer and inner SWNTs. 
Attaching a piezoelectric electrode to each edge, where 
only a single monolayer exists, we can measure the rela- 
tion between the interlayer motions and the total current 
forced to flow between the layers. The interlayer mo- 
tion and the interlayer force were investigated theoreti- 
cally using molecular dynamics^^ and experimentally by 
AFM and TEM.^^^ The interlayer conductances were 
also measured experimentally. ^^^^^ The relationship be- 
tween the interlayer conductance and the interlayer mo- 
tion can be used to construct nano electro-mechanical 
systems (NEMS), such as nano- mechanical switches^ 
and nano-displacement sensors. 

The interlayer bond between atoms a and /3 is repre- 
sented by the interlayer Hamiltonian element H^^fs- The 
interlayer motion influences the conductance through the 
change of these elements. Considering the interlayer 
Hamiltonian as a perturbation, we can show that the 
effective interlayer interaction is the product of the inter- 
layer Hamiltonian and the wave function. Though this 
effective interlayer interaction^^ was discussed in Ref.^^, 
its relation to the conductance was complicated. In the 
present paper, the perturbation formula simplifies this 
relation. The perturbation formula was not discussed in 
most of the preceding theoretical works about the con- 
ductances of TDWNTs^^^^. In Refs.2'^'2^ the perturba- 
tion formula was discussed but limited to the incommen- 
surate interlayer configuration. The present paper shows 
that the perturbation formula is effective even for com- 
mensurate TDWNTs, e.g. the (5,5)-(10,10) TDWNT, 
despite their larger interlayer interaction. In order to 
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FIG. 1: Telescoped double- walled nanotube (TDWNT). 



confirm the validity of the perturbation theory, the ef- 
fect of higher order terms , that was not discussed in 
Refs i^^i^^ , is also examined. 

Approximate analytical formulas in Refs.^^^^^ include 
fitting parameters other than the Hamiltonian. Though 
the fitting parameters are useful for precise reproduction 
of the exact results, they produce ambiguity about the re- 
lation between the conductance and the Hamiltonian. In 
the present paper, the perturbation formula is discussed 
for the interpretation of the relation between the Hamil- 
tonian and the conductance. Thus the fitting parameters 
are excluded. Since there is no fitting parameter, agree- 
ment between the perturbation formula and the exact 
results is limited. Nevertheless the perturbation formula 
is effective in the interpretation when it reproduces qual- 
itatively the relation between the Hamiltonian and the 
conductance. 

The paper is organized as follows. The tight-binding 
(TB) model used in the present paper is described in Sec. 
II. The perturbation formulas are derived in Sec. HI. The 
results of Ref.^^ are reproduced by the perturbation for- 
mulas in Sec. IV A. Corrections of the TB Hamiltonian 
suggested by Refs*^*^ are analyzed in Sec. IV B and IV 
C. Summary and discussion are shown in Sec. V. 



2 




12 3 4 



■>7 



->z 



a 



FIG. 2: Relationship between integer indexes (/, j) and coor- 
dinates iO^z). Triangles and circles correspond to odd / and 
even /, respectively. 



II. TIGHT-BINDING MODEL 

We considered the TDWNT composed of two armchair 
tubes, (no, no) and (n/, n/). The symbols O and / indi- 
cate the outer and inner tubes, respectively. As the in- 
terlayer distance must be close to that of graphite, only 
the case of no = + 5 was considered. Henceforth, 
the symbol /i indicates either O or /. The cylindrical 
coordinates (r, ^, z) of the atoms in the tube /i are 



RuL 



• 5^jA0, -j + 5^^o^z (1) 



with integers /, j, xij = 6/ — 5 — (—1)^+-^, the lattice 
constant a = x 0.142 nm, and i?^ = \/3anj^/{27r). 
Regarding the range of j, j > 1 for tube O and j < 2L 
for tube /, where L is the number of unit cells in the 
overlap region. As shown in Fig. 1, the overlap length 
equals (L — 0.5)a — Az. Figure 2 shows the relationship 
between the integers (/, j) and the coordinates (^, z). The 
geometric structure and the definition of (A^, A2:) are the 
same as in Ref.^^, where |A^| < 7r/no and \Az\ < a/4. 

The TT orbital at position ([!]) is denoted by and 
is assumed to be orthonormal, i.e., {ji' ' ^ j j) = 
'^1,1 '1 where 'bra' and 'kef are used to simplify 
the notation. The Hamiltonian H of the TDWNT is 
decomposed as = i^j^^ + Hq^ + V, where H^^^ and V 
correspond to the intralayer and interlayer Hamiltonian 
elements, respectively. The interlayer element between 
atom a = (O, j) and atom /3 = {I^l \j is represented 
by 



{W- 



■AW)e- 




cos Ua,/3 



■r <rc 
■r > Tc 



with 



inter- atomic distance 



(2) 



con- 



stants W 0.36 eV, d 0.334 nm, 0.045 nm, 

and cutoff radius Tc. The interlayer bonds were clas- 
sified as 'AA', 'BB', or 'AB' bonds. When either a 
or P had an interlayer bond shorter than ro with the 
third atom 7 (7^ Q^,/^) and r > ro, the bond be- 
tween a and P was classified as an AB hondi^^^ 
Here, ro = (a/V^) v'(aM)2TTT572^ = 0.34283 nm^ 
AV^ = -0.20 eV for AB bonds, and AW = for AA 
and BB bonds. The covalent bond character is the origin 
of the negative value of AW, the small cutoff radius rc, 
and the dependence on 0. The intralayer elements were 
—t = —2.75 eV between nearest neighbors, e for the di- 
agonal terms of Hq \ and zero otherwise. Even when no 
interlayer interaction exists, the linear dispersion lines of 
tube O shift from those of tube / due to the difference 
in curvature causing a — tt mixing. In the following dis- 
cussions, this shift is called 'intrinsic' and is distinct from 
the shift induced by interlayer interaction. ^^^^^^^^ The pa- 
rameter € was introduced to represent this intrinsic shift. 
The energv eI^I and the wave function |/i,(T, /c} of 

an isolated SWNT were obtained from H^^fi^cr^k} = 



-2t cos{ka/2) — ta -\- eS^^o 



l^^cr.k} = --^^exp(z/caj72)|/i,a,j)) 



and 



2n^ 



lM,Cr, j)) 



(3) 



(4) 



(5) 



1=1 



with the wave number k and the mirror symmetry a = 
±1. When 6 = and rc = (a/V^) ^(1.37)2 + (15/27r)2 = 

0.39085 nm,^ the total Hamiltonian H = Hf^ + V 

becomes equivalent to that of Ref^^. The reflection at 
the open edges (j = 1,2L) was neglected here, but will 
be considered in Sec. HI B. The ket |/i,cr, j )) defined by 
Eq. (jSj) win be used in Sec. HI B. 



III. DERIVATION OF THE PERTURBATION 
FORMULA 

A. Fermi's golden rule 

Since the interlayer Hamiltonian element in Eq. ([2]) 
is much smaller than the intralayer tt bonding t = 2.75 
eV, it can be considered as a perturbation. According to 
Fermi's golden rule (FGR), the probability of a transition 
caused by the perturbation V per unit time is 



P{0,(T,k\I,(j',k') = 



^\{0,cjMV\I.cj',k'}\^ 



(6) 
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The density of states with positive group velocity was 
derived from Eq. (|3]) as 



27r \dk 



= (^^/4t2-(^(M)+at-£(5^,o)2 ) (7) 



per unit cell of tube /i. Note that the wave function (j4]) 
was also normalized per unit cell. 

When the Fermi level is close to zero, the inter layer 
current / can be estimated to be 

a, a ' 

xP{0,a,k\I,(T',k') 



where 
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and 



j=ls=-l 



L-1 / (/) \ 2m 
m=0 \^<r / 



(15) 



(17) 



(18) 



where 



(9) 



Go = 2e^/h, di^^ = D^^^Ef), and V5 denotes the bias 
voltage. The Fermi wave number ki^"^ satisfies Ep = 
-2tcos{ki^^a/2) -ta ^ e, Ep = -2tcos{k''J] a/2) - ta ' 

and the positive group velocity, sin (/c^^^ a/2) > 0. As 
both the bias voltage V5 and the temperature were close 
to zero, the Fermi-Dirac distribution function difference 
fo ~ fi was replaced by eVbS{EF - E^^^) in Eq. ([8j). 

According to Landauer's formula, on the other hand, 
the conductance G = is determined by 



G — Go Terser 



(10) 



V;f, ^ {{0,a,3\V\I,a',3')) 
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{0,l,j\V\I,l',j'). (19) 



When j < or j ' > 2L + 1, Eq. i^9^ equals zero. 
Otherwise Eq. (fT9|) is determined by j — j ' and the 
parity of j. The cutoff distance Vc in Eq. (|2]) is so short 
that Eq. (fT9|) becomes zero when |j — j '| > 1. In Eq. 
(dH), {0,fcS^\cr|V|/,/c^^,\cr'} is resolved into the z-axis 
factor 5cr,cr ' and 6>-axis factor Acr^a ' with the boundary 
correction A^^^j] at z = (L — 0.5)a. The correction A^^^j\ 
however, is comparable to the ^-factor A^^a'^ while the 
z-factor 5cr,cr ' can become much larger than unity. Thus 
we can neglect A^^^j} in Eq. (p!l|) as 



where T(j^a ' denotes the interlayer transmission rate from 
|/, (t',/c'} to |0,cr, /c}. By comparing Eq. ([8j) to Eq. 
([To]), an approximate formula for the transmission rate 
can be obtained as follows. 
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T... ' = — I {O, a, 4""^ a ^ fe^',) } I (11) 



Here we concentrate our discussion into cases where \Ef\ 
and \e\ are much less than t, i.e., 



V3t 



(12) 



(13) 



andF^,,^/ -4/(3^2). 

Using Eqs. (Ill),® and (fTTjl . 



T , = 1- 
3^2 



4 _L 4 ,R 



(14) 



IB. 



(20) 



Because Bcr,-a <^ ^cr,cr5 ^+,- and are negligible 

compared to T+^+ and T_^_. Thus the following discus- 
sion will concentrate on the dominant transmission rates, 
T+,+ and T_,_. Since k^^^a - k^^^a 2e/{V^t), 
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When \e\L ^ >/3t , Eq. (|2T]) is approximated by 



3^2 



(21) 



(22) 



B. Green's function 

With the base set |/i,cr, j)) defined by Eq. ([5]), tubes 
O and / can be approximated by chains. The nonzero 
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FIG. 3: Schema of the reduced Hamiltonian when L — 2. 



intra-chain elements are ((/i, j| H^^^fi^j ± 1 )) = — t and 

{{li,j\H^^\li,j)) = ed^^o - crt. Here we suppress in- 
dex a to simplify the notation. As shown in Fig. [3] 
, he was introduced to cut away the artificial chains 
and form the open edges. The nonzero elements are 
((/,2L + l|/i,|/,2L)) = ((/,2L|/i,|/,2L + 1)) = t 
and {{O,0\hc\O,l))_= {{O,l\hc\O,0)) = t. The 
inter-chain element was defined by Eq. (p!9|) . 

The retarded Green's functions were defined with pos- 
itive infinitesimal 77 as ^ = {E ^ ir] — H^^^^ , g = 
{E^ir]-H^^^ -hcy\ G = {E^ir]-H^^^ -Vy\ 
and G = {E ^ if] - H^^^ - he - , where H^^^ = 



H^P + Hq^ . The inter-chain elements of g and g are where 



r(0) 



zero, while the intra-chain elements {{/J^^m \g \ii^n)) and 
{{^l,m\g |/i,n)) are denoted by (g^,)^^^ and {g^,)^^^, re- 
spectively. 

As was shown in Ref*^, 



I m—n I 



(23) 



where Eq. (p!3]) was used. Using Dyson's equation ^ 
5^ + 5^^c^7 we can derive 



Substituting m = 1 in Eq. (^^, 

is obtained. Using Eqs. (|23]), ([24l) and (|25]) . 

i9o)m,n = (9o)m,n " (5'o)^+n,0 ' 

The matrix gj is obtained in the same way as 

Dyson's equations G = g -\-g VG = g^^{V gY 
G = g^gVG = g Emi^dr indicate 



G 



and 



goVgi{l-'VgoVgi) 
goVgi {l-'VgoVgi)' 



(24) 
(25) 

(26) 

(27) 
and 

(28) 
(29) 



, respectively, where Gjj / = ((0,i | G |/ , j ' )) and 
Gjj / = ((O, j \G\I J' ))'. Using Fisher-Lee relation^^^^^, 
and Eq. (fT3]) . we can obtain the transmission rate T as 



(30) 



T = 3t^ G2L+i,o 
including the higher terms and 

T = 3t^ |G2L+l,o| 



(31) 



including both the higher terms and reflection at the open 
edges. The explicit relation of Eq. (j30j) and Eq. (|3T]) 
to the inter-chain elements Vjj > is summarized in the 
Appendix. In the Appendix, we can see that the first 
order term of Eq. (|3Q|) , 3t^| (^o'^5'/)2l+i I^' coincides 
with the FGR formula (dH). 



C. Expansion of Eq. ((30)) 

Equation (|3Q|) can be expanded as 



T = 3f 
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^=0 



and 



/5/ = —a I 




(33) 



(34) 



(35) 



(36) 



As the ranges of indexes of Eq. (|33]) are 1 < < 2L, 1 < 
j^' < 2L, -1 < < 1 and -1 < < 1,^^ the number 
of terms in Eq. (j33]) is 32^+^2^)^"^+^. Firstly we con- 
sider the case of £ = 0, i.e., k^^^ = k^^\ Among the 
22m+ij^2L)2^+^ terms, those satisfying the condition 



jo > Jo > ji > Ji > • • • > Jm-l > j^ 



(37) 



are dominant, because /3o + /5/ = 5^ + Y^^q^{^i + -^J) 
irrespective of indexes {juj'i) in these dominant terms. 
The other terms cancel each other because of their ran- 
dom phases Po ^ Pi- Since the number of the dominant 



5 



terms is 32^+i(2L)2^+V(2m + 1)!, Eq. (j33]) can be ap- 
proximated by 



(2m + 1)! 



A\A 



2m 



(38) 



where A was defined by Eq. (p!6|) . From Eqs. ([32]) and 
|), we can obtain 



(39) 



In order to discuss the case where A:^^) ^ k^^\ we refer 
to the dispersion relation of the DWNT. The DWNT can 
be approximated by the double chain of which the inter- 
chain Hamiltonian element V- ^ , is related to Eq. (p!9|) 
as V'^ y = {Vi^i^j'-j + F2,2+j '-j)/2. The dispersion 
relation of the double chain is obtained as^^^^^ 



■,/c ■ 2 



where 



2 1 

j = l s = -l 



(40) 



(41) 



Figure m shows that Eq. (|4Q|) coincides well with the exact 
dispersion lines. Comparing Eq. (|^T]) with Eq. (p!6]) . we 
can see that Aa- and Acr^a are essentially the same. In Eq. 
(|4Q|) . we can see that the intrinsic band shift s changes 
the total band shift - E'^'^ from \ A\ to Vp]^^^. 
Assuming the same effect of e on Eq. ([38|) . we can obtain 



2m+l 



-A(|^|2+e2)-. 



(2m +1)! 

From Eqs. and (021, we can also obtain 



■ sm 



(42) 



(43) 



When LV'IAP <^ |^| <^ |^|^ 

cides with Eq. i^. 



IV. ANALYSIS WITH THE PERTURBATION 
FORMULA 

In the following sections, the transmission rates at 
Ef = 0.1 eV 0.0364t are calculated from the condi- 
tioned transfer matrix (CTM)^^'^^ and the perturbation 
formulas (|20]) , (|30j) , (|3l]) with the common TB Hamilto- 
nian. In the CTM, the transmission rates were obtained 
from the 5* (scattering) matrix and the numerical errors 
were estimated to be ^.^^ KZ^j 5'*i5'j,/e) - ^i,k\, as the 
exact 5* matrix must be unitary. The estimated errors of 
the CTM in this paper were less than 4 x 10~^. Though 
the perturbation results were less accurate than the CTM 
results, they are useful in analyzing the CTM results. 
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FIG. 4: The dispersion relation of the DWNT with param- 
eters (no,n/, A6', Az) = (10, 5, -0.05a) and s (eV) =0, 
0.3. The triangles and circles show the exact dispersion rela- 
tion. The solid and dashed lines represent Eq. (|40p . 



A. Original TB model ( s = 0, rc = 0.39085 nm) 

Figures [5] and [6] represent + and T_ _, respec- 
tively, as a function of integer L. The overlap length 
(L — 0.5)a — Az was changed discretely and Az was fixed. 
The black solid lines, dotted lines, closed diamonds and 
open circles show Eqs. (|20]), dSO]), dH]), and CTM, re- 
spectively. In contrast to the monotonic increase of FGR 
formula ([20]) with L, CTM and Eq. §^ showed a rapid 
oscillation with a period close to 3a superimposed on a 
slower oscillation. Only a slower oscillation appeared in 
Eq. (|3Q]) because Eq. (|3Q|) does not include the refiec- 
tion at the open edge. To show the slower oscillation, the 
closed diamonds were connected with the dashed lines at 
intervals of 3a and the averaged CTM data defined as 
rj^iave) ^ ^rj.^_^ ^Tl^ ^l+i)/3 wcrc showu by the red 
solid lines. 

The effects of the structure parameters (L, no, Az, 
) were reproduced qualitatively by the first order for- 
mula (|2Q|) when Eq. (|2Q|) is less than unity, i.e., when 
L < V^t/\A\. The values of \A\/t are shown in Figure 
captions. Even when L > V3t/\A\, Eqs. (j30j) and (j3T]) . 
which include higher order terms, were effective, indicat- 
ing the validity of the perturbation formula. Note the 
scale of the vertical axis in Fig. [6l The transmission rate 
- of the TDWNTs becomes larger particularly when 
no = n/ + 5 = 10.^^ 



B. Modification of Tc 

When no = n/ + 5 = 10, the maximum conductance 
of the TDWNT is 2Go in the tight-binding (TB) calcu- 
lation, but only Go in the local density approximation 
(LDA) J^ii^ i21 The difference in the interlayer Hamilto- 
nian between TB and LDA is the most probable origin of 
this disagreement. To evaluate the order of the interlayer 
Hamiltonian elements in the LDA calculation, the ADF 



6 



CTM o averaged CTM 

Eq.(20) Eq.(30) 

Eq.(31) (a) (10,5) 0-0 




25 50 
(b) (10,5) 2 -(- 2) 





'/// / 

T/ 1 

iLi \ 


¥C\ / / 
■^'A\v / //^ - 


25 50 

(c) (14,9) 7-0 




6/ : 1 
"/r 1 


% 

\V /J 


//// 4. 





25 50 

the number of unit cells in the overlap region (L) 

FIG. 5: Inter layer transmission rates T+,+ between the + 
channels of {no^rii) ie-iz TDWNTs at Ef = 0.1 eV, where 
= 27Tie / (ISno) and Az = 0^^/40. The black solid lines, 
dotted lines, closed diamonds, open circles and red solid lines 
show Eqs. (Uni), (|30l), (|3l]), CTM and averaged CTM, respec- 
tively. The closed diamonds are connected with the dashed 
lines at intervals of 3a. The values of |A+,+ |/t are 0.1492, 
0.1435 and 0.1355 in (a), (b) and (c), respectively. 



calculation^^^^ with single zeta Is, 2s, and 2p orbitals 
was performed for a (10,10)-(5,5) DWNT composed of 
four unit cells. The structure is represented by Eq. ([T]), 
1 < j < 8, A6> = 0, and Az = 0.025a. Dangling bonds at 
j = 1,8 were terminated by hydrogen atoms with a bond 
length of 0.11 nm. In the {R^O^z) plane, the C-C-H angle 
is 27r/3 as is the C-C-C angle. Geometric optimization 
was omitted, since a slight change in the structure is not 
relevant to the order of the inter layer elements. With 
the TT orbital ifj^^ defined as ^/^t^ = cos^^2px + sin^^2py, 
Fig. [71 shows the interlayer ADF Fock matrix elements 
of the TT orbitals as a function of atomic distance. The 
TB elements used in Ref.^^ are also shown in Fig. [71 
for comparison. The intralayer elements between near- 
est neighbors were —4.73 ~ —4.49 eV in the ADF and 
-t = -2.75 eV in the TB. Thus, Fig. [3 shows that the 
interlayer elements normalized by the nearest neighbor 
elements were larger in the ADF than in the TB model. 
When we adjust the TB model to reproduce the LDA re- 
sults, the adjusted TB model needs to have contradictory 
features; larger interlayer elements and smaller interlayer 
transmission rates compared to the original TB. To re- 
solve the contradiction, we should notice that the ADF 
result showed no clear cutoff radius Tc in Fig. [71 Inspired 
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FIG. 6: The same calculations as in Fig. [5] except the trans- 
missions T-,- are between the — channels. In (c), Eq. ()20p 
and Eq. (|30p nearly coincide with each other. In the insets 
with the range 50 < L < 120, the data of Eq. ()3ip are omit- 
ted because they make the insets busy. The values of | |/t 
are 0.0167, 0.0262 and 0.0036 in (a), (b) and (c), respectively. 



> 0.7 

0.4 
^ 0.2 



0.3 



7 



interlayer atomic distance ( A ) 

FIG. 7: Interlayer Fock matrix elements between tt orbitals 
calculated by ADF. Interlayer tight-binding model based on 
Ref.^^ is shown as crosses for comparison. 



by these results, we discuss the effect of Tc in this section. 

Figure [51 shows T_ _ for = 0.385 nm, 0.39085 nm 
and 0.4 nm. We can see that T_^_ increases as Tc de- 
creases in the averaged CTM results (dotted lines). As 
this relation between Tc and T_^_ was reproduced by the 
perturbation formulas (sohd lines and dashed lines), it 
can be explained by the effective interlayer interaction 
- as follows. Interlayer bonds are drawn between 
atoms (0,/,l) and 2) in Fig. [9l when the corre- 




20 40 60 80 100 120 
the number of unit cells in the overlap region (L) 

FIG. 8: Inter layer transmission rate T_,_ of the TDWNT 
(no = m + 5 = 10, A6' = 47r/130, = 0) calculated at 
= 0.1 eV using Eq. (|20)) (solid lines), Eq. (|30)) (dashed 
lines) and the averaged CTM (dotted lines) for the shortened 
cutoff radius rc = 0.385 nm, the original one rc = 0.39085 nm 
and the lengthened one rc — 0.4 nm. The values of |A_,_|/t 
are 0.0271, 0.0160 and 0.0044 when rc = 0.385 nm, 0.39085 
nm and 0.4 nm, respectively. 



spending Hamiltonian elements (O, / , / 2 ) are fi- 

nite, or when the atomic distances are smaller than Vc 
The interlayer bonds with even or odd I ^ V are called 
even or odd bonds in the following discussion, and rep- 
resented by solid or dashed lines, respectively, in Fig. [9l 
where the parity of / and / ' is distinguished by triangles 
and circles, using the same representation of parity as in 
Fig. 2. 

The parameters ( no, A^,Az) were chosen to be the 
same in Figs. [9] (a) and (b) as in Fig. [51 Since T_^_ 
in Eq. (|2Q]) is reduced by different parities of / + / ' in 
Eq. (p!9|) . the reduction of T_,_ is determined by the bal- 
ance between odd and even bonds. This balance tends 
to be lost when the number of interlayer bonds per unit 
cell is small, as is illustrated by Fig. [9] (a); the even 
bonds are considerably longer than the odd bonds, al- 
though they are equal in number. On the contrary, this 
imbalance is redressed in Fig. [9] (b) by increase of the 
number of interlayer bonds caused by increasing Vc In 
cases where no > 10, on the other hand, the number of 
interlayer bonds is large enough to reduce T_^_ without 
increasing rc. This is illustrated in Fig. [9] (c), where 
no =ni^5 = 20, AO -7r/60, and Az = 0. Figure HO] 
shows |A_ _(rc = 0.4nm)p - |A_ _(rc = 0.39085nm)p 
for the (l6,10)-(5,5) TDWNTs of which the interlayer 
configurations are = ieir/Gb {ie = 0, 1, 2, • • • , 12) and 
Az = {i^ = 0,1,2,3,4). Figure ID] clearly in- 

dicates that increase of Vc tends to reduce lA-^-p and 
T_ _. 



(a) (^o-lO) 

=0.39085 nm 



(b) (^0=10) 

=0.4 nm 



j=l ♦-.^^^'^''^^^^^ 



(c) («o=20) =0.39085 nm 
0., 



J'=2 
7=1 



FIG. 9: Horizontal positions of triangles and circles indi- 
cate coordinates of atoms (0,/,l) and (/,/^,2) (a), (b) 
no = ni + 5 = 10, A6' 47r/130, 0, and (c) 

no = ni + 5 = 20, A6' -7r/60,Az 0. The inter- 
layer bonds were drawn when the inter-atomic distance was 
smaller than rc, so that the corresponding Hamiltonian ele- 
ments shown by Eq. ([2|) were finite. The cutoff radius rc 
was 0.39085 nm in (a) and (c), and 0.4 nm in (b). The inter- 
layer bond between atoms a and /3 was formed only when 

\0c. - Op\ < Oc = 2arcsin (^y^^An^r^ - X))/{12ninoa^) ). 

Here X = (75 + 7v^{j ' - jf)a^, j ' - j = 2 - 1 = 1, and 
Oc is shown by bars. 
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FIG. 10: Changes of \A-,-\^ of the (10,10)-(5,5) TDWNTs 
caused by increase of rc from 0.39085 nm to 0.4 nm. The in- 
terlayer configurations are A^ = ie'K /Q^ (i^ = 0, 1, 2, • • • , 12) 
and Az = i^a/40 (iz = 0, 1, 2, 3, 4). 



C. Modification of e 

Figure 6 of Ref.^^ indicates a close correlation between 
the intrinsic band shift \£\ and the suppression of the 
transmission rate T- -. These appeared in multi-band 
TB^^ and LDA, but not in single band TB Here we 
should distinguish \e\ from the total shift Y^|Ap + e'^ 



shown by Eq. (|^Q|) . Inspired by these results, we in- 
vestigate the effects of \e\ in this section. 

Figures [TTI and [T2l show the same calculations as Figs. [5] 

(b) and[6](b), respectively, except that the intrinsic band 
shift e was changed from zero to (a) 0.1 eV (b) 0.3 eV, or 

(c) 0.5 eV. The decrease in the CTM transmission rate 
with increasing e was reproduced qualitatively by Eq. 
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FIG. 11: The same calculations as in Fig. [5](b), except the 
intrinsic band shift e was (a) 0.1 eV, (b) 0.3 eV, or (c) 0.5 eV. 
The values of |yl+,+ |/t are shown in Table I. 



(|2Q|) . We can also see that precision of the perturbation 
formula was systematically improved by Eqs. (|3Q|) and 

m- 

Equation (|3Q|) almost coincides with Eq. (|43]) as is seen 
in Fig. [131 Thus Eq. (|43|) indicates that the first peak 
of Eq. ([30]) as a function of L appears at 



{L,T) 



(44) 



The first peak position of the averaged CTM is denoted 
by (^ctm, ^ctm) and compared to Eq. (|44|) in Table I for 
Figs. (b), il (b), H] and [H The peak height of Eq. 
(|44|) is lower than Tctm when 1^1 /|^| is large. Neverthe- 
less Eq. qualitatively reproduces the dependence of 
(Lctm, ^ctm) on £. When \£\ <C the intrinsic shift e 
exercises only slight influence over \A\. 



V. SUMMARY AND DISCUSSION 

Considering the interlayer Hamiltonian as a perturba- 
tion, we derived the first order formula ([20]) , the for- 
mula including the higher-order terms (|3Q|) , the formula 
including both the higher-order terms and reflection at 
the open edges ([3T]) . Expanding Eq. (|3Q|) . we can see 
that Eq. (|3Q|) is the essentially the same as Eq. (|^ . 
They were applied to TDWNTs composed of (no,no) 
and (no — 5, no — 5) armchair tubes. 

The perturbation formulas clarified the effects of the 
interlayer Hamiltonian on the interlayer conductance G. 




25 50 

(c) e = 0.5 eV 



o 

o 

,*4^'V--. 









25 50 

the number of unit cells in the overlap region (L) 



FIG. 12: The same calculations as in Fig. [6](b), except the 
intrinsic band shift e was (a) 0.1 eV, (b) 0.3 eV, or (c) 0.5 eV. 
The values of |A_,_|/t are shown in Table I. 



The product of the interlayer Hamiltonian and the wave 
function can be considered as the effective interlayer in- 
teraction because it determines the perturbation formu- 
las (|20l), (gal) and the dispersion relation ([40]). The 
effective interlayer interaction per unit cell and that 
per total overlapped region were denoted by Acr,cr ' and 
Acr^a 'Ba^a ' 1 respcctivcly. Here the parity a = +, — or 
a ' = +, — indicates whether the wave function of the 
outer or inner tube, respectively, changes its sign along 
the circumference. The first order (FGR) formula (|2Ql) 
was proportional to |Acr,cr ' Pl^cr,cr ' Because |5cr,-cr| 
was negligible compared to |5cr,cr|, we could neglect the 
inter-channel transmission rates T+^_ and 

The CTM transmission rate had a rapid oscillation su- 
perimposed on a slower oscillation as a function of L. 
Although Eqs. (|2Q|) and ([30]) could not reproduce the 
rapid oscillation, they approximated the long-period os- 
cillation. The rapid oscillation is reproduced by Eq. ([3T]) , 
which includes the reflection at the open edge. The first 
order formula ([20]) exceeds unity when < l^l, while 
the higher-order formula ([30]) never exceeds unity and re- 
produced qualitatively the first peak of the CTM results. 
The systematic improvement of the accuracy indicates 
the validity of the perturbation formulas. 

To represent the range L < Lmax where Eq. ([2Q|) 
reproduced the averaged C TM results, we define Li = 
t\/?>/\A\ and L2 = ntV^/ ^/4\A\^T4e^', Equation (|22|) 
reaches unity at L = Li and the first peak of Eq. (|3Q|) 
appears at L = L2. The upper limit of the effective 
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FIG. 13: Comparison between Eq. (|3Q|) (dashed lines) and Eq. 
(|43l) (solid lines) when (no, nj, AO, Az) = (10, 5, -Q.Q5a) 
and s (eV) =0, 0.1, 0.3, 0,5. We can see that Eq. gO]) almost 
coincides with Eq. (|43)) . The data of Eq. ([30]) are shown also 
in Figs. [5];b),[6i;b),[II]and[l2l 



g(eV) 



0.1 



0.3 



0.5 



0.1435 0.1435 0.1432 0.1430 



cr = +1 (Lctm) 



18.95 
(18) 



18.38 
(15) 



15.11 11.76 

(13) (14) 



{Tctm 



1.000 0.940 0.633 0.382 
(0.998) (0.807) (0.873) (0.570) 



\A-,-\/t 0.0262 0.0268 0.0279 0.0291 

, , 103.9 60.25 24.16 14.78 

2^|A_,_|2+£2 

-1 (Lctm) (90) (70) (28) (16) 



1.000 



0.351 0.062 0.025 



\A_ _|2+£^ 

(tctm) ( 0.918) (0.787) (0.149) (0.061) 



TABLE I: The first peak positions of the averaged CTM 
(Lctm, Tctm) in Figs. [5] (b), [6] (b), [11] and [12] are compared 
with Eq. (|44|) . The values of \A\/t are also shown. 



range, Lmaxi is classified according to |6:|/|74| as follows, 
(i) When > 1, Eq. ^ almost coincided with Eq. 

(|3Q|) and showed the underestimated peak heights. Nev- 
ertheless it reproduced well the period of the oscillation 
of the averaged CTM even when L > L2. In this sense, 
Lmax > L2. Figures [12] (b) and [12] (c) correspond to this 
case. {iVjLjnax = L2{< Li) when is comparable 

1.2). Figures [n](c) 

= L,{<L2), 



to unity but larger than y ^ — l(c 
and[T2](a) correspond to this case. {m)Lr, 



when \e\/\A\ < y x " 1- Figures [SE [3 [H] (a) and [H] 

(b) correspond to this case. 

Since the first principle results suggested the significant 
effects of the cutoff radius Tc and the intrinsic band shift e 
on the conductance, they were analyzed by the effective 
interlayer interaction Acr^aBa,a in Eq. (|2Q]) . The band 
shift l^l reduced the conductance because it lowered \B\. 
As Tc became longer, the number of nonzero terms in 
Eq. (p!9|) increased. When a = a' = — , the nonzero 
terms introduced by increasing Vc could cancel the terms 
already present in Eq. ([T9|) . Thus, a longer Tc could 
reduce and T_^_. 

In contrast to rigorous calculations, which involve com- 
plicated matrix inversion, the first order perturbation for- 
mula ([2Q|) involves a simple linear combination of the in- 
terlayer Hamiltonian elements. This enabled us to clarify 
the role of the interlayer Hamiltonian in the transmission 
rate. In addition to DWNTs, there are various other sys- 
tems composed of two monolayer subsystems; side-to-side 
contact of two SWNTs,^^!^*^^ an SWNT on graphene,'^^. 
and bilayer graphene.^^ By sliding one subsystem along 
the other, a generalized telescoped system can be ob- 
tained in which the interlayer bonds can be considered a 
perturbation. The perturbation formula is an important 
tool to analyze the NEMS formed by these telescoped 
systems. 
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Appendix 

Eqs. (|30|) and ^ are explicitly related to Eq. ([T9D 
in this Appendix. Only the case where a = a ' is consid- 
ered here and thus index a is suppressed in the following 
formulas as in Sec. Ill B. For example, the symbol t^^Jw 
is abbreviated as Vjj > . 



Eqs. (|3Q|) and (|3T]) can be represented by 



3r 



G2L+1,0 



and 



3t^ |G2L+l,o| 



2L 



2L 



(A.l) 



(A.2) 



with the following formulas 
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V3t(u;(0))"'''"'(5o^5/) 



2L+l,j' 



2L+1 

E 

io=i 



(A.3) 



y = (i + x(«)x(^))-\ 



(A.8) 



(A.4) 



^3,3' ^i,4L+2-i" 



(A.9) 



3,3 ' 



i{Vgi)3 



'^3,3+s 



S— — 1 ^ 



, (A.5) 



'2L+l,i' 



2L+1 

E 

i=i 



E2sm^jUjJ„ 



(A.7) 



and 



^3 ',3 - ^3 ',3 



3' -3 



(A.IO) 



For the inverse matrix calculation in Eqs. (jA.4p and 
dA^j, v^^^ and X(^) are considered as (2L + 1) x (2L + 1) 
matrixes in which indexes are restricted to 1 < j < 2L + 1 
and < < 2L. From Eqs. (US}, and 
we can obtain Eqs. (j32j), ()33p . (j34p . (j35j) and (j36j). 
Eqs. (jA.Sp and (|A.5p show that the first order term of 



Eq. (5^0^ ^1)2^+1^0 I 

the FGR formula ([14]). 
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